The shape and motion of a ruck in a rug 
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The motion of a ruck in a rug is used as an analogy to explain the role of dislocations in the 
deformation of crystalline solids. We take the analogy literally and study the shape and motion of 
a bump, wrinkle or ruck in a thin sheet in partial contact with a rough substrate in a gravitational 
field. Using a combination of experiments, scaling analysis and numerical solutions of the governing 
equations, we first quantify the static shape of a ruck on a horizontal plane. When the plane is 
inclined, the ruck becomes asymmetric and moves by rolling only when the the inclination of the 
plane reaches a critical angle. We find that the angle at which this first occurs is larger than the 
angle at which the ruck stops, i.e. static rolling friction is larger than dynamic rolling friction. 
Once the ruck is in motion, it travels at a constant speed proportional to the sine of the angle of 
inclination, a result that we rationalize in terms of a simple power balance. We conclude with a 
simple implication of our study for the onset of rolling motion at soft interfaces. 



Since Volterra's first mathematical description of dis- 
locations more than a century ago, the role of these lo- 
calized defects in the deformation of crystalline materials 
is well known pQ. To understand how these dislocations 
allow for slip along lattice planes, an oft-used analogy, 
attributed to Orowan pQ , connects the motion of an edge 
dislocation and that of a ruck in a rug. Just as it is easier 
to move a rug by having a wrinkle or ruck roll along it 
rather than by dragging the entire rug, it is easier for a 
crystal to deform by having a localized defect or dislo- 
cation glide or climb rather than sliding globally along 
an entire plane. When taken literally, the analog of a 
ruck in a rug appears in a variety of guises in scientific 
problems. For example, at interfaces between soft and 
stiff materials, e.g. rubber on glass, motion occurs by 
the generation and propagation of small wrinkles known 
as Schallamach waves [2], not by sliding. In thin films 
and filaments that interact with a substrate frictionally 
or adhesively, motion may again occur via the rolling of 
inchworm-like wrinkles [12 [4]. In flagellar axonemes, in- 
dividual microtubules often slide relative to each other 
while propagating microscopic rucks [5], while in cells, 
blebs, which are blisters where the cell membrane is de- 
tached from the underlying cortex can also propagate in 
a manner similar to that of rucks [6 j. Finally, earth- 
worms, leeches, and certain worms move using propagat- 
ing wrinkles that have been studied and mimicked [71 |8] , 
and surely played a role in creating the analogy in the 
first place. 

Motivated by these examples, here we study the statics 
and dynamics of a thin film or filament interacting with a 
substrate. In Figure [l] (a) we show the result of pushing 
a thin latex sheet (thickness h = 0.25 mm) lying on a 
rough plane till it buckles into a ruck that is stabilized 
by frictional forces. This "defect" may be characterized 
by the difference in the length between the sheet and the 
projection onto the substrate e, a quantity analogous to 
the Burgers vector for a dislocation (although it is not 
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FIG. 1: (a) A ruck in a rug on a flat substrate. The ruck 
height is A, while its length measured along the arc is S and 
its horizontal extent S— e, so that e is analogous to the Burgers 
vector. The latex rug is 0.25 mm thick, and the perspective 
view is a result of the tunnel view through the ruck, (b) The 
plane is at an angle of inclination, 0. The free-body-diagram 
in the upper right hand corner shows the local stress and 
moment resultants for a segment of length ds of the ruck, 
where <j) is the angle between the local ruck tangent and the 
plane. 



a discrete object here); indeed after the ruck has moved 
through the rug, the rug has been translated by a dis- 
tance e. When the plane on which the ruck sits is tilted, 
the ruck remains stationary until the inclination of the 
substrate reaches a critical value, and only then does the 
ruck move, and then does so at a constant velocity, shown 
schematically in Figure [I] (b) . These simple observations 
raise a number of questions associated with the shape 
and dynamics of a thin sheet in partial contact with a 
rough substrate, which we now address. 

Assuming the ruck to be inextensible, a good approxi- 
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FIG. 2: (a) The scaled ruck height, A as a function of the ex- 
cess length e, for different materials and thicknesses, h, (0.25 
mm (+), 0.5 mm (o), 0.75 mm (A), 1.0 mm (*), 1.5 mm 
(o)) and for the numerically calculated solutions (solid line) 
of ruck shapes collapses on to a single curve. All lengths 
are normalized by the elastic gravity length, £ g = (-j^)* • 
The inset shows the numerically calculated values (o) of A 
as a function of e on a log-log plot and compares well with 
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the scaling law A ~ given by (4). (b) Ruck shape as a 
function of the excess length, e; solid lines correspond to ex- 
perimental measurements, while the dashed lines correspond 
to numerical solutions of (l)-(3) 

mat ion when the latex sheet is compressively loaded, im- 
plies the kinematic relations between the instantaneous 
position of material points in the ruck x(s,i),y(s,i) and 
the orientation of the local tangent vector 0(s, t) rela- 
tive to the plane on which the ruck sits, parametrized in 
terms of the arc length coordinate s 

# )S = cos(/>, y yS = sine/). (1) 

Here dA/da = A iCL . Force balance parallel and perpen- 
dicular to the plane yields the equations of motion (see 
Figure [I] (b) for a pictorial derivation) : 

Fi,a + pghsinO + fi = 0, 
F 2)S ~ pgh cos + f 2 = 
M i8 - F 1 s'm(j) + F 2 cos<j) + m = 0. (2) 

Here F\ and F 2 are the integrated stress resultants in the 
film of density p, M is the torque resultant and g is the 
gravitational acceleration, and /i, f 2l m are the volumet- 
ric forces and torque at a cross-section (and include any 



potential contributions from inertia or other body forces) . 
Closure of the system of equations requires a specification 
of the torque M . Latex rubber is a viscoelastic material 
which we model as a simple Voigt solid, so that we can 
write M = EI a (j)^+nI a (j)^ s t, where I a is the area moment 
of inertia of the cross-section, E is the elastic modulus 
and fi is the viscosity of the material. Finding the shape 
of either the static or dynamic ruck requires the speci- 
fication of boundary conditions at the two contact lines 
that demarcate the locations where the ruck leaves and 
regains contact with the substrate, and read 

ar(0) = 2/(0) = 0(0) = >S (O) = y(S) = 0(S) = 0. (3) 

Here S is the length of the ruck that is free, and we have 
chosen the origin of the ruck to be the left contact line. 
The condition of the curvature vanishing at either end 
follows from the fact that there are no localized torques 
at the contact lines [TP] . 

For small amplitude rucks, the excess length e <C *S, 
and furthermore, the projected length of the ruck I ~ S. 
Then, geometry implies that the ruck height A = \/eZ, 
and its curvature n ~ A /I 2 . Comparing the elastic bend- 
ing energy, U e ~ Eln 2 !, with the gravitational poten- 
tial energy, U g ~ pghlA, together with the geometric 
relations k ~ A// 2 , / ~ h 3 yields the dimensional scal- 
ing laws I - e x l\Eh 2 lpgfl\ A - e 4 / 7 \Eh 2 / ' pg) 1 ' 7 , and 
the dimensional energy per unit length required to form 
a ruck scales as pghig 7 e 5 / 7 . Using the elastic gravity 
length, £ g = (^~)^ as a natural length scale in the prob- 
lem then allows us to write the dimensionless length and 
height of the ruck as 

l/£ g ~e$,A/£ g ~e$. (4) 

For rucks of large amplitude, we must solve the free 
boundary problem ([l][3| in the static limit (with fi = 
f 2 = m = 0), after normalizing all lengths by the elastic 
gravity length, £ g = stresses by N = pgh£ g , with 

the dimensionless variables x = xj£ g ,F\ = Fi/N,e = 
e/£ g , etc., and then drop the hats unless otherwise noted. 
In Figure [5] (a), we show the dimensionless height A as 
a function of the dimensionless excess length e, obtained 
by solving ([l][3| using a shooting method implemented in 
MATLAB. We also plot the experimentally determined 
height of the ruck for various values of e and see that the 
experimental and simulation data collapse onto a single 
curve. In the inset, we plot the same data on a log- log 
plot and see that our scaling law (El works for small to 
moderate values of e. In Figure T2] (b), we show the 
shape of the ruck for various values of e, determined ex- 
perimentally and numerically, and see that they compare 
well PZQ. 

When the plane on which the ruck is formed is tilted, 
the ruck becomes asymmetric, but does not move. Since 
only the gravitational energy scales differently, following 
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our earlier argument yields the scaling laws for the di- 
mensionless length of the ruck ljt g ~ e 1 / 7 (cos 6)~ 2 / 7 and 
its dimensionless height A/£ g ~ e 4 / 7 (cos#) -1 / 7 , while 
the dimensional energy required to create a ruck on a 
tilt scales as pghfj 7 e b / 7 (cos 6) A / 7 . To go beyond scaling 
arguments, we solve the equations (l)-(3) for various val- 
ues of the tilt and excess length e leading to the results 
shown in Figure [3] (a) along with the experimentally ob- 
served shapes, which compare well for moderate values 
of e but show substantial deviations for large values of e. 

When the angle of inclination is larger than a critical 
threshold, the ruck begins to move. It does so by rolling 
rather than sliding. Measuring e at the beginning of the 
run to that at its end, we find that it is conserved to 
within 2%. Starting with a stationary ruck in a rug on 
an incline, and increasing the angle of inclination slowly 
until it starts to move allows us to determine the criti- 
cal angle g at which rolling starts,. Once the ruck is in 
steady motion at a speed V, the incline is untilted slowly 
(relative to the movement of the ruck, i.e. 6 yt <C V/S) to 
determine the angle 9 S at which the ruck stops. In Fig- 
ure [3] (b), we see the marked difference between static 
and dynamic rolling friction characterized in terms of 
g and 6 S ; 6 g is an indicator of a critical torque which 
a ruck (of given shape and size) must overcome in or- 
der to begin rolling (a measure of static rolling friction), 
and S < g measures the dynamic rolling friction re- 
sistance, analogous to well known, but still incompletely 
understood difference between static and dynamic sliding 
friction. The existence of a critical inclination for the on- 
set of rolling is consistent with the existence of a critical 
torque due to microscopic interactions at the contact line. 
Interestingly, there is a similar threshold for the motion 
of dislocations pQ, but we will not pursue the possible 
mechanisms 

When 6 > 9 g , the ruck rolls down the inclined plane, 
quickly reaching a steady speed V and a steady shape 
that is similar to its static shape, at least for small ve- 
locities. In Figure [3] (c), we show the variation of V with 
the inclination for different values of e and observe that 
V ~ sin 6. Since the ruck moves by rolling rather than 
sliding, the gravitational power must be balanced by ei- 
ther air drag or dissipation within the ruck. For small 
amplitude rucks, balancing the viscoelastic power dissi- 
pation per unit length ph 3 (VA/l 3 ) 2 l [12 with gravita- 
tional power input per unit length phlVg sin#, we find 
that the dimensional speed V ~ pge~ 2 ^ 7 f^ 7 smO/ph 2 , 
i.e. smaller rucks travel faster than bigger rucks, with 
a speed proportional to the sine of the angle of the in- 
cline [13] . While the latter observation is consistent with 
our observations as shown in Fig. 3 (b), we do not see 
evidence that smaller rucks are faster than larger ones. 
We attribute this to the fact that experimentally, small 
rucks do not move until the angle of inclination is so large 
that a large part of the sheet slips before the ruck moves. 
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FIG. 3: (a) The ruck shape on an inclined plane with sin# = 

0. 3 shows the increasing asymmetry as e increases. Solid lines 
correspond to numerical solutions of (l)-(3), while the dashed 
lines correspond to experiments, (b) The ratio of the excess 
length to the contour length e/S vs. sin at the onset of 
motion with — 6 g and when motion stops with — S < 9 g , 

1. e. dynamic rolling friction is less than static rolling friction 
(c) Ruck speed V vs. sin# when > S for three values 
of e, 0.88,1.25 and 1.88, is consistent with the scaling law 
V ~ sin# determined by the balance between visco-elastic 
power dissipation in the latex and gravitational power input 
(see text). All experiments were done with a latex sheet of 
thickness h = 0.75 mm, and the shapes of the ruck (inset) 
correspond to 6 = 20°, 25°, 30°, e = 1.88. 



Thus, we are experimentally limited to studying a range 
of relatively large values of e when the shape of the ruck 
does not show the simple scaling behavior that led to this 
conclusion. 

To complete our investigation of the rolling ruck, we 
compare it to the rolling of a rigid wheel on a rigid plane, 
where the particle paths are simple cycloidal arcs sepa- 
rated by cusps. In Figure [4] (a), we show the path of a 
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FIG. 4: (a) The path of a point on the rug as a ruck moves 
it is plotted in white (*) on top of the time-lapse images of 
the ruck. In the inset is a plot of a point on the ruck as it 
approaches its zenith (*), compared the cuspidal curve y ~ 
x 2 ^ 3 that characterizes the kinematics of a point on a rigid 
wheel rolling on a rigid substrate, (b) For a slowly moving 
ruck, the dynamic shape of a ruck is close to its static shape 
(dashed line) along with the particle path of a point on it 
(inset, solid line). The experimental curves correspond to 
e = 1.25, = 20°, V = 0.7 m/s while the numerical curves 
correspond to e = 3.5, = 17°. 



particle in the rug that is transported by the ruck. In 
the inset of Figure [4] (a) we plot the local variation of the 
trajectory of a material point in the ruck in the neigh- 
borhood of the cusp, consistent with the curve r] ~ £ 2//3 
that characterizes the cuspidal profile for a particle on 
a cycloidal trajectory. In Figure [4] (b), for comparison, 
we show the particle path using a numerically calculated 
shape of the stationary ruck. We see that while parti- 
cle paths in a rolling ruck are similar to those in a rolling 
wheel in the vicinity of three points (when it starts, stops 
and near its zenith), in general the rolling wheel and 
rolling ruck are quite different, given the deformable na- 
ture of the latter. 

Our study has answered some basic questions about 
the shape and motion of a soft ruck in a rug, but 
leaves us with the question of understanding the tran- 
sition from statics to dynamics in the ruck. Here, we 



wish to simply point out that this requires us to ac- 
count forthe inclusion of the effect of frictional inter- 
actions at the elastic contact lines. Then the contact 
line torque EI(j) jS (0) ^ EIcj) jS (S) ^ in general, so 
that global torque balance on the ruck then implies 
M r = EI(/) i8 (0) - EI^ S (S) + pgSlp(e/S,6) + 0, where 
p(.,.) is a dimensionless function that characterizes the 
shape of the ruck, and the inclination of the plane. At 
the onset of rolling, M r /pghSl = M*(0, e/S, ...) = M*, 
a scaled threshold torque. Determining M* in terms of 
the various interactions at the contact line constitutes a 
basic question in the dynamics of soft rolling and indeed 
in the context of all the applications mentioned in the 
introduction. However, just as phenomenologically de- 
rived velocity weakening laws allow us to make progress 
on certain sliding friction problems by circumventing the 
question of the onset and stoppage of sliding at interfaces, 
for rolling at interfaces one might be also able to do the 
same, given that we now have experimental evidence for 
the existence of a similar scenario. 

Acknowledgment As this work, begun six years ago 
in Cambridge, was being completed, we became aware of 
very recent work addressing the inertial dynamics of the 
motion of a ruck by D. Vella, A. Boudaoud and M. Adda 
Bedia, who we thank for their preprint. 



* llm@seas.har vard.edu 
[1] F. Nabarro, Theory of dislocations, Dover, 1982. 
[2] A. Schallamach, Wear, 17, 301-312, 1971. 
[3] J.H. Gittus, Phil. Mag., 31, 317-329, 1975. 
[4] K. Kendall, Nature, 261, 35-36, 1976. 
[5] S. Aoyama and R. Kamiya, Biophys. J., 89, 3261-3268, 
2005. 

[6] G. Charras, T. Mitchison and L. Mahadevan, Biophys. 

J., 94, 1836-1853, 2008. 
[7] E.R. Trueman, Locomotion of soft bodied animals, 

Arnold, 1975. 

[8] L. Mahadevan, S. Daniel and M. Chaudhury, PNAS, 101, 
23-26, 2003. 

[9] Landau, L. Lifshitz, E., Theory of Elasticity, Pergamon 
Press, New York, p. 52, 1970. 
[10] The condition 4>, S (S) = is automatically satisfied as 
can be verified by multiplying the last equation in Q by 
<p 8 , and integrating the result in light of the other two 
conditions. 

[11] There is a critical value of the excess length e c above 
which the ruck loses symmetry and tilts to one side, a 
question that we will not address further here. 

[12] Here we have used the relation that the power dissipated 
is ph(JiK t ) 2 l, with K t ~ Vk s ~ VA/l 3 . 

[13] For small amplitude fast moving rucks, where the resis- 
tance arises from air drag, balancing the gravitational 
power phlVg sin 6 with the power to overcome air drag 
p f V 3 A 2 yields V ~ (pgh sin 0/ p f ) 1/2 (£ g /e) 3/14 . 



